#! /usr/bin/env python

import pyx
import pyfeyn
from pyfeyn.user import *

fd = FeynDiagram()
processOptions()

in1 = Point(-4,  2)
in2 = Point(-4, -2)
out1 = Point(4, -2)
out2 = Point(4,  2)

in_vtx = Vertex(-2, 0, mark=CIRCLE)
out_vtx = Vertex(2, 0, mark=CIRCLE)

fa1 = Fermion(in1, in_vtx).addArrow().addLabel(r"\Pq")
fa2 = Fermion(in_vtx, in2).addArrow().addLabel(r"\Paq")
bos1 = Photon(in_vtx, out_vtx).addLabel(r"\PZ")
bos2 =  Photon(out_vtx,out2).addLabel(r"\PZ")
higgs=Higgs(out_vtx,out1).addLabel(r"\PHiggs")

fd.draw("WZ_Higgs_production_Feynman.pdf")

